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AUTOMATIC CONTROL SYSTEIVIS SATISFYING CERTAIN GENERAL CRITERIONS ON 

TRANSIENT BEHAVIOR ' 

By Aaron S. Boksenbom and Richard Hood 


SUMMARY 

An analytic method for the des-ign of automatic controls is 
dtceloped that starts from certain arbitrary criterions on the 
behavior of the controlled system and gives those physically 
realizable equations that the control system can follow in order 
to realize this behavior. The criterions used are developed in 
the form of certain time integrals. 

General results are shown for systems of second order and of 
any number of degrees of freedom. Detailed examples for 
several cases in the control of a turbojet engine are presented. 

» 

INTRODUCTIOiN 

In the past several years, there has beeu increasing develop- 
ment and interest in automatic control; in the fields of gun 
direction, guided-missile control, and control of gas-turbine 
('iigines, for instance, where very refined and accurate 
controlled behavior is required, need still exists for further 
development of the methods of controls analysis and design. 
' Recent developments in this field have been concerned 
chiefly with the problem of control analysis both in the 
realm of linear systems (reference 1) and in the realm of 
nonlinear systems (reference 2). These analytical works 
answer questions concerning the behavior of a given system 
or the effect on its behavior of changing certain constants 
in the system. 

Another problem of equal and, in some cases, greater 
importance is that of control synthesis. Work on this prob- 
lem seeks to find the method of design when certain cri- 
tcrions concerning the behavior of a controlled system are 
given. The method decided upon should give aU aspects 
of the system; for instance, where the system should be 
linear, what the general configuration should be, and what 
the precise values of all the constants should be. 

This STOthesis problem has hardly been broached in liter- 
ature. The use of analysis as a design procedure offers a 
[lartial solution to this problem in that the analysis of a 
large number of cases may reveal, coincidently, one that 
satisfies the desired criterions of controlled behavior. Such 
a method is, at best, long and tedious and almost always 
would result in compromises because the systems chosen to 
be analyzed would probably be such that they could never 
satisfy all the desired criterions. 

A method for designing a linear system to satisfy certain 


special criterions when operating on a random input is 
developed in reference 3. This method is applicable as an 
addition to a control system whenever random external 
disturbances are involved. A partial solution to the sjm- 
thesis problem is developed in reference 4 in satisfying the 
criterions of noninteraction for systems with many degrees 
of freedom. 

An analysis made at the NACA Lewis laboratory during 
the latter part of 1950 and presented herein __develops a 
rational method of control synthesis that starts from any 
arbitrary but physically realizable criterions and results in 
the equations for the best system that satisfies these cri- 
terions. As is shown, the nature of the criterions, in general, 
requires minimizing certain time integi'als by using the 
calculus of variations and the methods .developed are an 
application of the calculus of variations to the problem of 
control synthesis. 

A careful scrutiny is first made of the whole problem, fol- 
lowed b}- a development of general results. These general 
results are then applied in examples to the design of turbojet- 
engine control systems. In general, the methods used vary 
according to the order of the differential equation describing 
the plant, the process, or the system being controlled and 
according to the number of degrees of freedom or independent 
variables. Detailed analyses are presented for application 
to a system of first order and of one degree of freedom. 
General results for systems of second order and any number 
of degrees of freedom are developed in the appendix. 

SYMBOLS 

The following symbols are used in this report: 

a,b,c gas-tuihine-engine characteristic constants 

C constant 

E function of X 

F, H functions of y and y 

'F(ti) indicates F[y(U), yit)] 

Fs,{y,y) indicates partial differentiation with respect 

to y 

f,fr arbitrarj" function 

fufz, . . . weighting functions used for gas-turbine con- 
trol criterions 

G, Gi, Gz functions of y 

g function used for gas-turbine-engine surge 

criterion 


’Supersedes NACA TN 2378, “Automatic Control Systems Satisfying Certain General Criterions on Transient Behavior” by Aaron S. Boksenbom and Richard Hood, 1951. 
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Subscripts: 

0 

f 

Superscripts: 

* 


function used for gas-turbino-eiigine blow- 
out criterion 
controller gain 

temperature limits for gas-turbine engine 
deviations of gas-turbine-engine speed, com- 
pressor-discharge pressure, characteristic 
temperature, and fuel flow, respectively, 
from values at some common equilibrium 
condition 

actual engine. acceleration 

exponent 

time 

time at end of transient 

variation in time at end of transient 

independent variables, functions of time 
small number 
arbitrary constants 

gas-turbine-engine time constants for re- 
sponse to temperature 
transient time constant of controlled system 

initial condition of variable 

final value of variable 

setting or desired value of variable 

indicates case diiferent from optimum 


The dot indicates dilTerentiation with respect to time. 

The prime indicates difi'erentiation v'ith respect to the argU' 
nient shown. 


ANALYSIS, 


SURVEY OF PROBLEM 


Control problem.— An important aspect of the control- 
synthesis problem is a clear definition of the criterions of 
desired controlled behavior. If a variable y is to be con- 
trolled, a reasonable criterion is that the time integral of 
some function of y is to be a minimum or a constant; that is, 

/(y) dt=constant or minimum (1) 

Ju 

or 

(y—ys)^df= constant or minimum (2) 

Equation (2), for instance, weights the error in y as the square 
and according to the time duration of that eiTor. . Another 
type of criterion may be that which requires a certain time 
duration to be a minimum or a constant; that is,' 



df=constant or minimum 


( 3 ) 


F. The, tise of a single criterion, such as equation (1), will 
'usually yield /(y)= constant. This result is reasonable be- 
cause jiy) can usually be made identically a constant if no 
additional criterions are imposed on other variables in the 


system. Certain limiting conditions usually exist on other 
variables in the system and these conditions must be included 
in the origmal criterions. 

Thus, a possible criterion could be written 


for 



—ysYdt = minimum 


rh 

/(s) constant 
Jo J 


( 4 ) 


If, for instance, y=engine speed and 2 = characteristic teim- 
perature of a gas-turbine engine, the criterion of equation 
(4) states That it is deshed to design a control system such 
that, for a particular value of a temperature integral, the in- 
tegral of the speed-error squared is a minimum. This cri- 
terion may be used if, for instance, it is laiown that an over- 
temperature condition can be tolerated for a certain period 
of time and it is desired to keep the average speed error at 
a minimum during a transien t. 

The general theory will show tliat as many criterions as 
desired of the type shown in equations (1) to (4) can be in- 
cluded together and a control system can be derived that 
automatically satisfies all these criterions simultaneously. 

Another aspect of the control criterions is the end con- 
ditions of the integrals of equations (1) to (4). The time in- 
terval for which th(!se integrals are to be a minimum or a 
consttyit must be chosen. A reasonable time interval is any 
duration during which essential external disturbances are 
constant and during which the system to be controlled moves 
from one. essential level of operation to another. The essen- 
tial external disturbances are those that cannot be imme- 
diately corrected by the control system. If an essential 
external disturbance were allowed within the time interval 
of the criterions, no physically realizable system could be 
expected to anticipate this disturbance so as to behave prop- 
erl}’' before this disturbance occurs. An essential level of 
operation is any specific condition of only those variables 
that must be continuous. It will be shown that the essential 
level of operation appears as a natural boundary condition 
for the type of criterion used. In the case of a turbojet 
engine the transient behavior of which can be described by 
a first-order differential equation, tlie engine speed deter- 
mines the level of operation. If a lag exists in the fuel 
system or. between temperature and engine speed, then both 
engine speed and acceleration are required to describe the 
essential operating level of the engine. 

Analytic problem,— The control system resulting from any 
design method must be physically realizable. There are 
two aspects to this problem. Fust, it is possible to set down 
criterions that are not realizable with any system or are 
incompatible with each other. If such criterions are used, 
the unrealizability will appear either as a requirement on 
the control to look ahead into the future or as an inability 
to satisfy the boundary conditions of some differential equa- 
tion. In most cases, a clear understanding of the criterions 
used and of the system to be controlled will indicate incom- 
patibilities of this sort. 
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The second aspect of physical realizability is purely mathe- 
matical. It is desired to derive a description (a differential 
equation.) of the control or the controlled system that satis- 
fies the criterions of control and all the necessary boundary 
conditions that arise in the derivation of this equation. 
Although the mathematical solution of the problem may be 
any derivative or integral of this differential equation, the 
physical solution of the problem requires the differential 
equation that itself satisfies the boundary conditions and for 
which no undetermined constants of integration exist. Thus, 
such forms as 

y=Cx ) 

and > (5) 

U=(yx ) 

are not necessarily interchangeable as descriptions of some 
part of a controlled system because the forms differ by an 
undetei'mined constant of integration. For stable linear 
systems, the effect of this constant becomes vanishingly 
small; for the general nonlinear systems presented herein, 
however, this constant must be considered. 

Stability problem. — The requirement of stability is a 
special criterion that does not enter into the main body of 
the methods of this report. It may enter in the final steps of 
the method where tlie final differential equation describing 
the controlled system may be the integral of a higher-order 
differential equation that satisfies the necessaiy'- boundary 
condition for stability. In addition, it is always necessary 
to add to the controlled system a stability de\rice that does 
not affect the behavior of the system as far as satisf 3 dng the 
other criterions is concerned. This device can be described 
as follows: 

y=Q when y=y, (6) 

or, for a second-order system, 

y=0 when y=y, and y=0 (7) 

GENERAL THEORY AND RESLTLTS 

It has been shown that the criterions for control can be 
developed in the foUowing forms: 


\f{y)dt 



and so forth. If, for such a list of criterions,. one of the 
integrals is to be a minimum under the condition that the 
other integrals are to be constant, it is sufficient (reference 5) 
to make 


r ^f{y)dt+\i r \y — y,ydt+ 

Jo Jo 

*dt=minimum (9) 
or 

*L/(i/) + >^i(l/ — ( 2 ) -hX 3 ]df= minimum (10) 

The X’s are arbitrary' constants that enter into the control 
system as the adjustable parameters and are precisely de- 
termined by the choice of values that the constant integrals 
are to have. 

The technique of the X multipliers is widely used for prob- 
lems of this ty'pe where one condition is to be a minimum 
under other restrictive conditions. Indeed, the conditions 
need not be in integral form and any’- functional or differential 
relation among variables can be handled in a similar manner 
(reference 6). 

Equation (10) can be made very general when all possible 
restrictive conditions are included. In the final equations, 
which are derived later, if any' one criterion is not to be used, 
then the corresponding X— >0. If any of the criterions is to 
be zero (as for the case of a variable having an absolute limitl,_ 
then the corresponding X— ^ ■= . 

If, for this development, the system to be controlled is of 
first order and of one degree of freedom (has one independent 
variable), then the variables y and 2 are such that 2 = 2 ( 1 /, y). 
Equation (10) can then be written, in general, as 


f ^ F(y,7j)dt=mm\mum (11) 

Jo 

In equation (11), Fis a continuous function of y and y, and 
■y is a continuous function of time. 

The calculus of variations (reference 5) is used to determine 
y as a function of time such that the integral of equation (11) 
is a minimum; that is, if the solid curve of figure 1 makes 
the integral a minimum, any other curve (such as the dashed 
one) win make the integral equal to or greater than the first 



Figubs 1 niustratioii of curves for minimization of 
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integral. If the two curves are verj^ close, the condition 
that the two integrals are equal will make the integral of 
equation (11) stationary (maximum, minimum, or inflection 
{7oint). 

The problem -is similar to that of finding the maximum or 
minimum point on a curve by setting the first derivative 
equal to zero. WTiether a maximum or minimum point 
exists is decided lyv the second derivative at that point. In 
the variational problem, proving a true minimum involves 
taking the “second variation.” As in the problem of finding 
a minimum on a curve, the second variation proves, at most, 
a local minimum and not an absolute minimum. For the 
specific examples discussed herein, an absolute minimum is 
proved by another method. In manj’ cases, the pli3'sical 
meaning of the equations will indicate the existence of a 
unique minimum that obviates going further than the 
“first variation.” 

On figure 1, the curve for minimizing the integral of equa- 
tion (11) is the solid curve having the value y{t) at an\' time. 
An\^ other curve differing from it b\^ a small amount is shown 
as the dashed curve having the value t/(i) + e hy{t) at anj' 
time, where e is a small number and hy{t) is an arbitrary' con- 
tinuous function of time. The condition for the integral of 
equation (11) to be stationaqv is 

F{y-\-t&y,y^i5y)dt=0 ( 12 ) 

at 

€ = 0 

where ^ [5?/(f)] = ^y{t) (reference 5). The time duration for the 

int(fgral of equation (12) is such as to start at some definite 
time (f=0) but not to end at a definite time but rather 
along some curve y=J{t) (fig. 1), in order to allow the proper 
boundaqv conditions of moving from one essential level of 
operation to another, as previouslj- discussed. 

Performing the operation indicated b}" equation (12) leads 
to 

1^"' |^Vf52/df+F(fi)5fi=0 (13) 

Integrating the second term of equation (13) by parts gives 

Jo ^ ^~~dt -fF(h)5fi=0 (14) 

Because by is an arbitrarj^ function, the integral and the 
boundary-condition terms must vanish separately. The 
geometry' at the end condition {t=ti, fig. 1) gives 

^y{U)=[J'{U)-y{U)\ bh (15) 

as €— »0. The t wo conditions that follow from equation (14) 
tile 11 become 

(16) 

and 


The time interval during which the criterion of equation 
(11) is to hold is considered as that during which the system 
moves from one essential operating level to another; in this 
case, from one definite value of y to another definite value 
of y. Thus, 

5i/(0) = 0) 

and > (18) 

/'(h)=o) 

Equation (16) is satisfied onl^- if the integrand is zero, and, 
because 0, the two conditions of equations (16) and 
(17) become 

(19a) 

and 

F{U)=Fy{U)y{ti) (19b) 

Equation (19a) need not hold at f=0 because 6y(0) = 0. 
The onh' condition that need hold at t=Q is that Ej(0) is 
finite, which will be true if y is continuous at f — 0. At 
the start of a new transient, y, F, Fy, and Fy maj' be dis- 
continuous, whereas at other points (0Ot<h) Fy will be 
continuous because of equation (19a). 

Equation (19a) is the differential equation for the y{t) 
that satisfies the original criterion of equation (11). The 
plu’sical answer to the problem is the first integral of equa- 
tion (19a), which satisfies the boundary-condition equation 
(19b). This solution is 

F{y, y)=yFy{y, y) ( 20 ) 

and whenever y, Fj, and so forth are continuous, 

y\Fy-J^{Fy)^=Q ( 21 ) 

b\' differentiating equation (20) with respect to time. Thus, 
either y—0 (which is true onl\' during stabilipv) or equat ion 
(19a) is satisfied. 

Thus, equation (20) is the description of that ph\'sicall\' 
realizable s\^stem the behavior of which wiU automaticalh’' 
and simultaneous!}' satisf}' those criterions included in the 
function F during that time interval for which the external 
disturbances are constant and during which the sj'stem goes 
from one operating level to aiqv other operating level. A 
stability device must be added to the s}'sLem; the description 
of such an ideal device is 

?/=0 when y=y, (22) 

An additional condition must be met if y is discontinuous 
in the interval In this case, Fy must be continuous 

during the discontinuity in y\Fy, however, will usuall}' be 
discontinuous when y is discontinuous. This discontinuit}' 
in Fy usuall}' means that some essential external disturbance 
has entered the s_vstem and the point of discontinuity must 
be the start of a new time interval. . _ 


bt,{F{U)+FM Lf(f0-i/(h)]}-Ei(O)52/(O) = O (17) 
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APPLICATION TO TURBOJET-ENGINE CONTROLS 

In the usual case of designing turbojet-engine controls, 
tlie engine speed that sets the essential operating level of 
the engine and, in the main, other pertinent characteristics 
such as thrust, is to be set or controlled. Limiting conditions 
of tlie engine are those of_ overspeed, overtemperature, 
compressor surge, and rich burner blow-out. The following 
criterions on the behavior of this engine are typical: 

rh 

Ji{N—Ni)df for speed control 

.h 


f 2 (AO dt 

Jn 

S,{T)dt 

Jo 


for speed overshoot 

for temperature overshoot 
and undershoot 


fi[P—g{N')\dt for compressor surge 


[P — h (AOl dt for blow-out 


r 

Jo 


dt 


for rise time 


(23) 


Figures 2(a) to 2(e) show the nature of the functions appear- 
ing in expressions (23). The variable involved is essentially" 
weighted by the function shown and integrated with respect 
to time. The quantity P—g(N) is the amount the 
compressor-discharge pressure exceeds the safe pressure for 
surge and g(N) (shown in fig. 2(f)) is the compressor-discharge 
pressure for each engine speed at a safe value below sui^e. 
Rich burner blow-out can be handled in a similar manner 
(k(N) is shown in fig. 2(g)). The rise time is the total time 
for the system to move from one essential operating level 
to the other. 

The linearized engine characteristics can be e.xpressed, by 
assuming first-order behavior, as foUows: 

T=aN'-ra(rN (24a) 

P=bN+cT (24b) 

Thus, the integrals of expressions (23) are of the form 


r 

./o 




(25) 


where/ is a continuous function of A^and N, and N (barring 
impulsive temperatures and the like) is a continuous function 
of time. 

Speed control; case A. — If only the error m speed control 
is considered important, the criterion becomes 

//A'’— A/) minimum (26) 


X 


where F=fi(N—N,). Equation (20) becomes 

/(A— A/)=0' ■ 

and from tlie nature of,/ (fig. 2(a)), / 

A^=A/. 


( 27 ) 




/.frj 



(c) 




(a) Sp«ed rontrol. 

(b) Speed OTershiwt. 

(0) Temperature OTerahoot and undershoot. 

(d) Compressor surge. 

Ce) Blow-out. 

(f) Compressor line below surge, 

(g) Compressor line below blow<iut. 

Fioijee 2.— .Arbitrary weighting lune.tions for various control criterions and pertinent engine 
characteristics involved. 


This result means that, in the absence of other criterions on 
the engine behavior, this speed control should keep speed 
error identically zero, which is physically realizable only in 
the sense of allowing infinite temperatures and the like to 
keep the speed error identically zero. This result, howe\'er, 
is inconsistent with the previous development of equation 
(20) in that A'^is now a discontinuous function of time and the 
time interval of the integral of equation (26) is zero. This 
instance is actually' a trivial case of the general problem. 
The result (equations (27)) does indicate that a criterion like 
that of equation (26) must be accompanied by an additional 
criterion (temperature overshoot, for instance) to give a 
physically realizable sy'stem. 

Speed control with temperature-limiting criterions; case 
B. — If the error in speed control and the overshoots and 
undershoots in temperature are to be considered the primary 
criterions of control, then from equations (10) and (11), 


£ ‘ L/iUV-A^) + dt = mi 


minimum 


(28) 


where F=fi{N'—Ns) + \fs(T). Equation (20) then becomes 
/ (zV- A/) + \UT)=\aaNMT) (29) 
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and file ideal stability device is such that 

iV=0 when N=Ng (30) 

Equations (29) and (30) describe the. complete control sys- 
tem. 

In figure 2(c), it is convenient to let/ 3 (T) = (T— X 2 )" for 
r>Z, and MT) = (T-Li)” for T<U In general, the 
power n should be ^ 1 , because, when n<l, T may be infinite 
and of such nature as to make iV discontinuous and physically 

('ll 

unreal even though jz{T)dt is finite. In the examples nf 

this report, n=2 Mid fi{N— Ns) = {N—N,Y . 

A variety of methods of settmg up the control system to 
realize equation (29) exists. IMien the preceding e.xpressions 
for'/i and /3 are used, equation (29) can be put in the con- 
venient form 


A 

where, for acceleration, when N<^Ns, 

iV>0andX=X2 

and, for deceleration, when iV>iVs, 
•Y<0 and X=Xi 


adNe 



When N < N„ N > 0 and L-Lz 
When N ^ N^, N < 0 and L 'Z, 


Fir.tiEE 3 .— Schematic diagram of control for turbolet engine for ease of speed control with 
temperature-limiting criterions (case B) . Stability device must be added to this system. 

Fa 

A schematic diagram of the system is shown in figure 3, 
where equation (31) is considered to give a desired N for anj" 
value of N, Ns, and X. Consistent with equation (31), a 
right-triangle construction is used to give a desired N. The 
actual Ns can be made verj" close to N by using a high-gain 
proportional controller, as shown. Provision must be made 
to change the sign of N and the value of X w'hen N—Ng 
changes sign. In addition, the stability condition requires 
a provision for making iV =.0 whenever speed error is very 
small or zero. 

The use of a high-gain proportional controller, which fol- 
lows from the requirement that N ma}^ be discontinuous, 
means that the fuel-flow rate required may be Infinite if lags 
exist in the fuel sj^stem or in the feedbacks. But, as no 
criterion has been set on fuel-flow rate, this requirement does 
not violate the original criterions. If necessary, however, a 
criterion on fuel -flow rate may then be added to equation 
(28). Even though a criterion on temperature is being satis- 
fied, no direct measurement of temperature is used in figure 3 . 



Actually, the equation for temperature (equation (24a)) is 
used as an indirect indication of temperature. 

The control sj'stem of figure 3 has one adjustable param- 
eter X. For any value of X, this system will, for the value 
of integral temperature-overshoot squared obtained, give the 
minimum value of integral speed-error squared. The value 
of X determines the actual value of the integral temperature- 
overshoot squared. 

The integral temperature-ovez'shoot squared as a func- 
tion of X is shown in figure 4 for the special case where 
aN,=L; that is, acceleration or deceleration to the speed 
that coiTesponds to limiting temperature. In this case, the 
system of figure 3 becomes a simple first-order lag system 
and equation (31) becomes 

E{L — aN) — atjN (32) 

where 



Figifre 4.— Various control parameters as functions of E for spued control.with temperature' 
limiting criterions (case B) when accelerating to limiting temperature. 


In figure 4, the iiitegi’al speed-error squared, the maximum 
temperature, and the time constant for this transient are 
also shown as functions of X. The curve showhig the mini- 
mum rise time for the corresponding temperature integi'al 
is included for a comparison that will be discussed laler. 
The equations for the curves of figure 4 arc 
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(time constaiit) 



(33) 


'Hie left sides of these equations have been put in dimension- 
less form. The maximum temperature T^ai occurs at the 
hemnning of the transient. The time constant r is that for 
the controlled system and is shown compared with the engine- 
time constant rr. 

From figure 4, the value of X is chosen as a compromise 
between the various quantities of equations (33). For£'=l 
(\=oDjj the temperature does not overshoot, the speed 
integral is 0.5, and r=o-. As E increases (X decreasing), 
the temperature integral and the maximum temperature 


a<SN 



Fh'rrRE 5.— Phase plane showing dviiamies of controlled system for speed control with 
temperature-limiting criterions (case B) where g^=I. 
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increase, whereas the speed integral and the time constant 
decrease. A value of (or a’X= 1) appears to be a good 

compromise and is used for the subsequent discussion. 

The behavior of the system of figure 3 can be best seen by 
drawing the trajectories in the phase plane shown in figure, 
5 where acN is plotted against aN for lines of constant 
aN, according to equation (31). Lines of constant temper- 
ature are 45° parallel lines in this plot and the lines of T=L 
are shown. Each trajector\' intersects, and is tangent to, 
the line T=L at N=N,. Figure 5 completely describes the 
transient behavior of the system. For any starting point 
anywhere on the plane (for instance, point A), the s}’stem 
will automatically move the operating point to that tra- 
jectory corresponding to the lYj that exists (point B) and 
then along this trajectory to the point T=L (iV=iYj, point 

C) , and finally the stability condition will enter to move the 
operating point along the solid vertical line to N=0 (point 

D) . 

The time sense for these transients is obtained by solving 
the differential equation (equation (31)) for the speed and 
the temperature time responses. The equations for these 
solutions are as follows; 


cosh E 


(34) 



Figcee 6. — ^Typical transient of controlled system for speed control with temperature- 
limiting criterions (case B) where 


These transients are shown in the general case in figure 6 for 
a step increase in Ay. Maximum temperature occui's at 
f=0 and the temperature overshoot decreases as A'' increases, 
and when N=N„ T=L and t=0. The stability condition 
then causes T to drop to its equilibrium value at this point. 
The time scale shown in figure 6 corresponds to the specific 
relative values of the ordinates shown. 
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Minimum rise time and temperature-limiting criterions; 
case C. — In order to obtain a minimum rise time for a con- 
stant temperature integral, the requirement is that 

[1 d-X/sCT)] minimum (35) 

where F=l-fX/ 3 (r). Equation (20) becomes 

l+X/3(T) = Xa<TiV//(r) (3G) 

and the ideal stability device is such that 

iV=0 when N=Ng (37) 

where Ng is the desired stable N. Equations (30) and (37) 
describe the complete control sjj'stem. Equation (36), and 
therefore the control system proper, does not include and, 
in this case, the stability device is independent of the control 
system proper. 

When the same/ 3 (T) is used as was used in the previous 
case, equation (36) can be put in the convenient form 

i+(i-aA0'=aVW2' 

A 

where, for acceleration, when A’’<A''j, 
iV>0 and 

and, for deceleration, when N^N„ 

N<C0 and L—L^ 




Figure 8,— Various control parameters as function of X for minimum rise tima with toinpora- 
ture-limitiiig criterions (ease C) when accelerating to limiting temperature. 


tion or deceleration to the speed that corresponds to limiting 
temperature. In addition, the rise time, the ma.ximnrn 
temperature, and the integral speed-error squared arc shown 
for this transient. A curve showing the minimum speed 
integral for the corresponding temperature integral (from case 
B) is also shown and will be discussed later. The equations 
for the curves of figure 8 are as follows: 


P'kjure 7— Schematic diagram of control for turbojet engine for mmimuin rise time with 
temperature-limiting criterions (ease C), StRbility device must be added to this sytem. 


A schematic diagram of this system satisfying equation 
(38) is shown in figure 7. This control system is the same as 


the previous system (fig. 3) except that in figure 7, 

A A 

replaces ^ — p-A on one leg of the right-triangle construction 


and the stability condition must be imposed outside the 
control system, as A(, does not enter into the criterion of 
control (equation (35)). 

The system of figure 7 has one adjustable parameter X, 
which, as before, sots the precise values of the integrals 
entering into the criterion (equation (35)), as well as all other 
behavior characteristics. 

The temperature integral is shown in figure 8 as a function 
of X for the special case where aNf=L; that is, for accelera- 


1 Y 

V Jo \.L gA'q/ 




{L—aNo)^^ 
•if dt=sinh"‘ [(L— aA^o) yx] 

O’ Jo 


T —T 

■*- max 


L-aN 


0 y'X {L — aN^ 


ro/ A^-A"^ V 

<rjo \L—aN^ 


dt _ vh+X(L — «A^o)^ 
2{L— aA^o)y'X 


1 


; sinh — a 


Y (39) 


2X(L-aA'o)' 

Again X is to he chosen as a compromise between these 


various quantities. In this case, -p=.has the units of 
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temperature and (-5— a-Vg) is the pertinent dimensionless 
[tarameter. For a fixed X, tiie initial speed lYo v'ill therefore 
cliiuige the time integrals and an increased L—aNo 'wfill 
decrease the temperature integral and the maximum tempera- 
ture and increase the rise time and the speed integral. 


acfN 



Fioore 9.— Phase plane showing dynamics of controlled system for minimom rise time with 
temperature-limiting criterions (case C) where (i— oA'o) = ^8. 


The dynamic behavior of the system of figure 7 is shown in 
the phase plane in figure 9. This figure is a plot of aaN 
against aN, according to equation (38), for various values of 
X. Lines of constant temperatures are 45° parallel lines and 
the lines T=L are shown. For any fixed X, only two trajec- 
tories are followed, one for accelerating and one for decelerat- 
ing. From any starting point on the plane, the system will 
automatically move to that trajectory corresponding to ac- 
celeration or deceleration and will move along this trajector 3 ^ 
until the stabHitj' condition enters (at N=N,) to make iV=0. 

The dependence of the time integrals on iYo maj" require X 
to vary with Na- In figure 9, A'o’s corresponding to each X 
are shown such as to keep X(L— axYo)^=8. The value of 
this parameter was so taken as to have the temperature 
integral in this case equal to that of the pre ious case for 
purposes of comparison. 


The time sense for these transients is obtained b\' solving 
the differential equation (38) for the N—t and T—t tran- 
sients. The equations for these solutions are 


/X (uiY.— L)=sinh ( — \-C 


^\(T-L)=e G+c) 


(40) 


These transients are shown for the general case in figure 10 
for a step increase in N,. The temperature will jump to 
some value above L, but in this case, unlike the previous one, 
the temperature continues .to increase as the speed increases. 
Alaximum temperature occurs at the end of the transient. 
Whenever N=N, (two such conditions are shown), the 
stabilit}' condition takes over and T is reduced to its equi- 
librium value. The time scale shown on figure 10 cotve- 
sponds to the specific relative values of the ordinates as shown. 



F ir.CEE 10.— Typical transients for controlled system tor minimum rise time with temperature- 
limiting criterions (case C). 


Comparison of cases B and C. — In ease B, a s^'stem was 
devised that, for a constant value of the integral temperature- 
overshoot squared, gives a minimum integral speed-error 
squared. In case C, a st’stem was devised that, for a con- 
stant temperature integral, gives a minimum rise time. For 
case B, the rise time is hot a minimum but can be compared 
with the minimum rise time of case G. For the special 
transient of accelerating to T=L, the s^'stem of case B 
reduced to a first-order s^’stem, the time constant of which 
is shown in figure 4. Because five time constants are con- 
sidered as the rise time of an exponential, the time constant 
of case B should be compared with one-fifth the minimum 
rise time of case G. In figure 4, for corresponding values of 
the temperature integral, one-fifth the minimum rise time 
from case G is plotted. Figure 4^ now shows that the rise 
time for case B is about twice the minimum possible rise time. 

At corresponding values of the temperature integral,. 



N—NfY dt for caseG can be compared with the minimum 


possible value of this quantity’ from case B. The minimum 


value of 


J* (N — NfY dt is plotted in figm-e 8. 


It is seen 
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tliafc, for the same temperature integral, (N—N;Y dt when 

minimum rise time is obtained is about 120 percent of the 

minimum possible value of j' (N—NfY dt. 

Proof of absolute minimum . — As previously noted, the 
minimization of equation (11) involves not only the condition 
that the first variation be stationary, which leads to equation 
(20), but also the condition that the second variation be 
positive. This second condition, however, would prove at 
most a local minimum for the condition of equation (11). 
A special proof of an absolute minimum is shown as follows: 
If equation' (11) is written in the form 

f ’^i^(i/, i/)^=minimum (41) 

Jiio y 

and y is considered as a function of y for the minimum con- 
dition and y* is considered as a function of y for any other 
possible case, then the condition for an absolute minimum is 
that 



or 

('‘\F{y,r)-tF{y,y)W,>Q (43) 

ji'o L y jy 

where dyly* is a positive differential. If the function F is 
considered in a general form, quadratic in y, as 

F{y,y) = Q{y)-^yGi{y)+y^Q.,{y) ' (44) 

then equation (43) becomes 

(45) 

The equation of the control system (equation (20)) becomes, 
for the form of F assumed in equation (44), 

fG.2{y)^G{y) (46) 

Using this expression for G{y) changes equation (45) to the 
form 

r G,{y)(.r-yy%>0 (47) 

Ji'o y 

For cases B and C, G 2 =\d‘a^ where X> 0. Thus an absolute 
minimum is proved for those two cases. 

Degree of minimum; case B.— The left side of equation (47) 
can be used as a measure of the deviation from optimum 
conditions when the methods of this report are not used. If 
the two cases {N for the optiihum case and N* for any other 
ease) are compared at the same value of the temperature 
integral, the left side of equation (47) is the difference between 
the integral speed-error squared for any case and the mini- 


mum possible value of this quantity. The ratio of this 
deviation to the minimum value becomes 

Fract.onal increase p,/p-TV dv 

m speed mte«ral= (jVr^) 7T ‘ ' 

for the transient of acceleration or deceleration to limiting 
temperature. The coefficient of the integral of equation (48) 
for the value of X previously chosen (a^X=l) is 2^l2. 

Case B can also be considered as giving a minimum value 
of the temperature integral for any definite value of flie speed 
integral. If two cases (iV for the optimum and N* for any 
other case) are compared at the same value of the speed 
integral, the left side of equation (45) is proportional (o the 
difference between the integral temperature-overshoot 
squared for any case and the minimum possible value of this 
quantity. This deviation can be written 

Fractional increase in 2E p/Z T* — T V ^1/ cqq') 
temperature integral \T/— To/ <^y* 

for the transient of acceleration or deceleration to limiting 
temperature. The coefficient of the integral of equation 
(49) for the value of X,previously chosen (a^X=l) is 16.4. 

Degree of minimum; case C. — If the two cases iVand N* 
are compared at the same value of the temperature integral, 
the left side of equation (47) is the difference between the 
rise time for an}" case and its minimum possible vmlue. This 
deviation can be written 

Fractional increase X(l— oiVo)^ pi/ T*— T V dy 
in rise time ~sinh“‘.^;'X (i — aiVo)Ji'o \Tf—TJ ay* 

(50) 

for the transient of acceleration or deceleration todimiting 
temperature. The coefficient of the integral of equation (50) 
for the value of .\/x(L— oA^o)=2>/2 previously chosen is 4.54. 

Case C can also be considered as giving a minimum value 
of the temperature integral for any definite value of rise time. 
If the two cases N and N* are now" compared at the same 
value of rise time, the left side of equation (47) is proportional 
to the difference between the integral temperature-overshoot 
squared for any case and the minimum possible value of this 
quantity. This deviation can be written 

Fractional increase m Jvo\ig—J- ;/ ay* 

temperature integral ~.^/x(Z;-aAro) 2+ 1 

(61) 

for the transient of acceleration or deceleration to limiting 
temperature. The coefficient of the integral of equation (51) 
is 16.4 for the value of ^[\{L—aNd) previously chosen. 

GENERAL SUMMARY 

I\Tien the criteiions on the behavior of a controlled system 
can be expressed in certain general forms, as 
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( 8 ) 


where the time interval is taken as any duration during which 
essential external disturbances are constant and during which 
the system moves from one essential operating point to 
another, the optimum control can be considered as one that 
minimizes one of the integrals of equations (8) while main- 
taining the other integrals at prescribed values. The 
analytical problem, according to the calculus of variations, 
reduces to the following equation: 

( \f(y)dt + \i r \y—y,ydt+ 

Jo Jo 

k" f '/o( 2 )dt-rX 3 f Vf=minimum (9) 
Jo ” Jo 

For general first-order systems, where z=z(y,y), equation 
(9) reduces to 

f’h 

•^(y, y)df=minimum (11) 


where F=f{y) + \i{y—y,Y-\-\J^{z)-\-\z. The equation nec- 
essary for the control s^-stem to satisfy equation (11) and 
all the boundary conditions becomes 

F(y,y)=yFj{y,y) (20) 

This equation should be followed by the control system 
proper. In addition, a stability device must be added to the 
system, the idealized characteristics of which would make 
y = 0 when y=y,. 

The arbitrary multipliers X are then found by evaluating 
the integral criterions involved in F. The transient behavior 
of the system derived is. found by solving the differential 
equation (20). The degree of the minimum or the amount 
suffered when any other control system is used was evaluated 
for the special cases considered. A summary of these devel- 
opments follows for a special form of the F function, quad- 
ratic in y, where 

F (y ,y) = Q (y) + y (? i(y) + r ’(?2 (y) (4 4) 


Control-system equation. — For F in the form of equation 
(44), the control-system equation (20) becomes 


Giy) 

G-iiy) 


(52) 


The function GJy) does not affect the control system and 


the control-system equation gives an explicit expression for 
y as a function of y. 

Evaluation of integrals. — If an integral of a function 
H(y,y) is any one of the criterions to be considered, it can be 
evaluated as follows: 


P‘i7(y,y)df= 

Jo Jro 


A^)J 


V G(y) 


(53) 


Thus the integrals can be evaluated without solving the 
differential equation (52). 

Transient behavior.-— The differential equation (52) can 
be easilv’ solved as follows: 



dx=t—to 


(54) 


Degree of minimum. — Equation (47) was derived to prove 
the absolute minim um for the special cases considered and 
can be used to evaluate the degree of the minimum found. 
Thus, 


rmy,r)dt*- rmy,y)dt^ G,(y)[r(.y)-m?J^ 

Jo Jo jsfo y ly) 

(55) 


where H(y,y)dt 


is to be a minimum or maximum. 


If ^2(y)>0 for yo^y^y/, then an absolute minimum is 
obtained; if (? 2 (y )<0 for yo<y<y/, an absolute ma.ximiini is 
obtained. Equation (55) also indicates that the degree of 
minimum (or ma.ximum) varies with the magnitude of G^iy). 
The y* in the denominator indicates that for any deviation 
|y*— yj in y, it is better to err on the side of a larger |y*|. 


SUMMARY OF RESULTS 

A rational analytic method- for the design of automatic 
control systems was derived. Criterions on the behavior of 
the controlled system were developed in the form of certain 
time integrals. When any of these arbitrary' but physically 
realizable criterions were used as a starting point, those 
equations that the control system must follow were derived. 
The criterions developed required the minimization of certain 
time integrals using the calculus of variations. The method 
gave not only a description of the behavior of the controlled 
s 3 ’stem but also those ph\-sically realizable equations that the 
control system can follow in order to realize this behavior. 
General results were shown for systems of second order and 
of anj’ number of degrees of freedom. 


Lewis Flight Peoplulsion L.abor.atory 
N.ation.al Advisory Committee for Aeronautics 
Cleveland, Ohio, October 11, 1950 
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APPENDIX 

SECOND-ORDER SYSTEMS OF SEVERAL DEGREES OF FREEDOM 


It is beyond the scope of this report to detail the cases of 
higher-order systems and those having several independent 
variables. The general equations for these cases will be 
developed and it is expected that subsequent reports will 
cover their applications in detail. 

For the case of a second-order system with two degrees of 
freedom, equation (11) becomes 


2 )df=mini 


minimum 


(Al) 


where y and z are independent functions of time. The 
condition to satisfy equation (Al) is 


d 

del 


at 


F{yFtby,yFeiy,yFe^y,zFehz, 

z-\- e&z,z-\- ehz)dt=Q 
(A2) 

e = 0 


The time duration of the integral of equation (Al) is that 
beginning at a definite time (i=0) but not ending at any 
definite thne but rather along some curves: 
g=g^(t), and z=gi{t) (fig. 11). The functions 5y and hz are 
arbitrary and independent functions of time. 






Figure 11.— Illustration of curves for minimization of | Fiyjy.VyZ.z.z} dt. 

JO 


Performing the operation indicated by equation (A2) 
gives 

r {FyBy-\-Fiihy-{-Fij^y^Fi5z-TFi5zFF-ib^dt-{-F{t-^hti=Q 

(A3) 

After integration by parts, 


dt 


(A.)-t- 


dt 


~2 (Pi) iz dt^(^ibyFF^5y—^^ Fydy'^ ^ -)-F(ti)St 


Fi&z+F-iSz 




= 0 


Jo 


1+ 


(A4) 


As before, the integrands of the integrals and the boundary- 
condition terms must vanish separately. From the geometry 
at the end condition fig. 11), 

5y(ti) = \Ji(ii)-yiii)] 


(A5) 


SyiU)=[f2(ti)S(ti)] sh 
^z(U)=[gi(ti)~2(ii)] 5^1 

6i(ti) = [^/(fi) — z(ti)] 5ti j 

The three conditions from equation (A4) then become 

[F-yFiFy^(F,)-yF,-zF,+z^^(F^)-z^^^^^ 

fi%) (A,-) ^ (A6) 

g/(tO [f, (Fs)]^^^ +^/(f-i) F-i(Q + 

SyiO) (Af)]^^+«l/(0) A’^(0) + 

52(0) (Fs)2^^-f 5i(0) Fm = 0 

The fii'st two of equations (AG) are the differential equa- 
tions that satisfy the original criterions of equation (Al). 
The phj'sical solution to the problem is the pair of solutions 
of these equations that satisfies the boundaiy-condition 
equation of equations (A6). If the fu’st of equations (A6) 
is multiplied by y and the second, by z and the equations are 
added, an exact derivative is formed, the integral of which 
is 

F-yFy+y ^ {Fi)-yFy-zFi+z ^ {Fi)-zF-i=C (A7) 
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(AlO) 


A study of the boundarj'-condition equation of equations 
(A6) then shows that the pliysically reasonable boundarj- 
conditions should be as follows: 

If ov FyT^K, then 5i/(0) = 0 and /i'(h) = 0 

If then 5y(0)=0 and 

If and Fi^K, then 5z(0) = 0 and^i'(fi) = 0 

If FzT^O, then 5i(0) = 0 and .f 2 '(ii) = 0 

Then in equation (A7), 6’=0 and the final solution to the 
problem of equation (Al) for the boundary conditions of 
equations (AS) is 

F-ijFi+y (Fi)-yFy-zFi+z^^(F,)~m=Q' 

and either 

f\~(Fy)-{-^AFs) = 0 >- (A9) 

or 

The meaiiincr of the boundary-condition equations (AS) 
is to define the original criterions for that duration during 
which the system moves from one essential operating level 
to another. Thus, if all conditions of equations (AS) must 
hold, the system goes from one definite y, y, z, and z to any 
other definite y, y, z, and i. The first difl'erential equation 
of equations (A9) would be of third order and the second or 
third equations of tA9) would be of fourth order. If equa- 
tions fA9) are integrated to obtain a pair of second-order 
differential equations having three constants of integration, 
the choice of these tkree constants can then determine a 
desired end point, that is, the values of y, z, and i at some 
final y. 

I'he physical solution to the problem, then, is the pair of 
secoiid-onler differential equations that are solutions of 
equations (A9j with the constants of integration evaluated 
so that the system goes through some desired end point. 
This end point must be such as to allow the possibility of 
stability. Such an end point may be written as follows: 


y=0, z = z„ and i=0 when y=y, 

which gives three conditions for the evaluation of the tkree 
constants of integration. A stability device must still be 
added to the system so that, at the point when equations 
(AlO) hold, 

y=o ) 

and V (All) 

z=0 ) 

Equations (A9) are symmetric in the variables y and z, 
which indicates the nature of the extension for more inde- 
pendent variables. Thus, if a third independent function 
w exists, the original criterion would be written 

j A(y,y, i/, z,z, A to, li?, fp) dt=minimum (Al2) 

This condition is satisfied under boundary conditions similar 
to equations (A8), where two additional conditions are added 
on the variable w, and the following equations describe the 
controlled system: 

F—yFj -h y ^ (Ff) - j/Ef — lEid-i ^ (Fj) —zFi— 
wFw-tVO ^ {Fii)—wFw=0 
and any two of the foUowing three equations: 

F.-^(Fi)-fg(F£) = 0 

The three equations of equations (A13) can then be mtegrated 
to give three second-order differential equations where the 
five constants of integration are evaluated so that the system 
goes through the desired value of y, z, z, w, w for some final 
value of y. 




REFERENCES 


1. .Anon: Theory of Servomechanisms, Hubert M. James, Nathaniel 

B. Nichols, and Ralph S. Phillips, eds. McGraw-Hill Book Co., 
Inc., 1947. 

2. Minorsky, N.: Introduction to Nonlinear Mechanics. J. T\'. Ed- 

wards, 1947. 

3. Weiner, Norbert: Extrapolation, Interpolation, and Smoothing of 

Stationary Series. The Tech. Press (M. I. T.) and John Wiley 
& Sous, Inc., 1949. 


4. Boksenbom, Aaron S., and Hood, Richard: General Algebraic 

Method Applied to Control Analysis of Complex Engine Types. 
NACA Rep. 980, 1950. (Formerly NACA TN 1908.) 

5. Bolza, Oskar: Lectures on the Calculus of Variations. G. E. 

Stechert & Co. (New A'ork), 1946. 

6. Lanczos, Cornelius: The Variational Principles of Mechanics. Univ. 

Toronto Press, 1946. 



